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We perform calculations of the dependence of nuclear magnetic moments on quark masses and 
obtain limits on the variation of (mq/ A.qcd) from recent measurements of hydrogen hyperfine (21 
cm) and molecular rotational transitions in quasar absorption systems, atomic clock experiments 
with hyperfine transitions in H, Rb, Cs, Yb^, Hg^ and optical transition in Hg^. Experiments with 
Cd"*", deuterium/hydrogen, molecular SFe and Zeeman transitions in '^He/Xe are also discussed. 
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I. INTRODUCTION 

Interest in the temporal and spatial variation of major 
constants of physics has been recently revived by astro- 
nomical data which seem to suggest a variation of the 
electromagnetic constant a = /tic at the 10~^ level for 
the time scale 10 biUion years, see (a discussion of 
other limits can be found in the review Q and references 
therein). However, an independent experimental confir- 
mation is needed. 

The hypothetical unification of all interactions implies 
that variation of the electromagnetic interaction constant 
a should be accompanied by the variation of masses and 
the strong interaction constant. Specific predictions need 
a model. For example, the grand unification model dis- 
cussed in Ref. predicts that the quantum chromody- 
namic (QCD) scale Aqcd (defined as the position of the 
Landau pole in the logarithm for the running strong cou- 
pling constant) is modified as follows: SAqcd/^qcd ~ 
34(5a/a. The variation of quark and electron masses in 
this model is given by Sm/m ~ 70 Sa/ a. This gives an 
estimate for the variation of the dimensionless ratio 



S{m/AQCD) ^ ^^5a 

{m/KqcD) Oi 



(1) 



This result is strongly model-dependent (for example, 
the coeflticient may be an order of magnitude smaller 
and even of opposite sign 0). However, the large co- 
efficients in these expressions are generic for grand uni- 
fication models, in which modifications come from high 
energy scales: they appear because the running strong 
coupling constant and Higgs constants (related to mass) 
run faster than a. This means that if these models are 
correct the variation of masses and the strong interaction 
scale may be easier to detect than the variation of a. 

One can only measure the variation of dimensionless 
quantities and therefore we want to extract from the 
measurements the variation of the dimensionless ratio 



rriq/ Kqcd - where ruq is the quark mass (with the depen- 
dence on the renormalization point removed) . A number 
of limits on the variation of rUq/ Kqcd have been ob- 
tained recently from consideration of Big Bang Nucle- 
osynthesis, quasar absorption spectra and the Oklo nat- 
ural nuclear reactor, which was active about 1.8 billion 
years ago (see also |i|lSlii|ll|l3). Below 

we consider the limits which follow from quasar absorp- 
tion radio spectra and laboratory atomic clock compar- 
isons. Laboratory limits with a time base of the order 
one year are especially sensitive to oscillatory variations 
of fundamental constants. A number of relevant mea- 
surements have been performed already and even larger 
numbers have been started or are planned. The increase 
in precision is happening very fast. 

It has been pointed out by Karshenboim 14] that mea- 
surements of ratios of hyperfine structure intervals in dif- 
ferent atoms are sensitive to any variation of nuclear mag- 
netic moments. First rough estimates of the dependence 
of nuclear magnetic moments on rriq/ Kqcd and limits 
on the variation of this ratio with time were obtained in 
Ref. j5|]- Using H, Cs and Hg+ measurements we 
obtained a limit on the variation of rriq/ Kqcd of about 
5 • 10"^'^ per year. Below we calculate the dependence of 
nuclear magnetic moments on rriq/ Kqcd and obtain the 
limits from recent atomic clock experiments with hyper- 
fine transitions in H, Rb, Cs,Yb"'",Hg+ and the optical 
transition in Hg+. It is convenient to assume that the 
strong interaction scale, Kqcd, does not vary, so we will 
speak about the variation of masses (this means that we 
measure masses in units of Kqcd)- We shall restore the 
explicit appearance of Kqcd in the final answers. 

The hyperfine structure constant can be presented in 
the following form 



A = const X 



i{Za)][^i- 



(2) 
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The factor in the first bracket is an atomic unit of en- 
ergy. The second "electromagnetic" bracket determines 
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the dependence on a. An approximate expression for 
the relativistic correction factor (Casimir factor) for an 
s-wave electron is the following 

Prel — J-. 5 7T , (3) 
7(472 - 1) 



where 7 = i/l — [Zaf- and Z is the nuclear charge. Vari- 
ation of a leads to the following variation of F^^i tlSj : 



K 



= A 

Frel a 



(4) 



(5) 



72(472 - 1) ■ 

More accurate numerical many-body calculations ^3 of 
the dependence of the hyperfine structure on a have 
shown that the coefhcient K is slightly larger than that 
given by this formula. For Cs (Z=55) K— 0.83 (instead 
of 0.74), for Rb ii:=0.34 (instead of 0.29) and finally for 
Hg+ K=2.2S (instead of 2.18). 

The last bracket in Eq. Q contains the dimcnsionless 
nuclear magnetic moment (i.e., the nuclear magnetic 
moment M = M 2m c ) ' electron mass rrie and proton mass 
•nip. We may also include a small correction arising from 
the finite nuclear size. However, its contribution is in- 
significant. 

Recent experiments measured the time dependence of 
the ratios of the hyperfine structure intervals of i^^Hg+ 
and H [3, ^^^Cs and ^^Rb [H and the ratio of the 
optical frequency in Hg+ to the hyperfine frequency of 
^^■^Cs In the ratio of two hyperfine structure con- 

stants for different atoms time dependence may appear 
from the ratio of the factors Frei (depending on a) as 
well as from the ratio of nuclear magnetic moments (de- 
pending on TTLq/ Kqcd)- Magnetic moments in a single- 
particle approximation (one unpaired nucleon) are: 



ioi j = 1 + 1/2. 



2(J + 1) 



{-gs + (2j + 3)gi) 



(6) 



(7) 



for j = I — 1/2. Here the orbital g- factors are gi = I for 
a valence proton and gi = for a valence neutron. The 
present values of the spin g-factors, gs, are gp = 5.586 for 
proton and (?„ = —3.826 for neutron. They depend on 
rUq/ Aqcd- The light quark masses are only about 1% 
of the nucleon mass (m^ = (m„ + md)/2 « 5 MeV) and 
the nucleon magnetic moment remains finite in the chi- 
ral limit, rriu = md = 0. Therefore, one might think that 
the corrections to gs arising from the finite quark masses 
would be very small. However, through the mechanism 
of spontaneous chiral symmetry breaking, which leads to 
contributions to hadron properties from Goldstone boson 
loops, one may expect some enhancement of the effect of 
quark masses [Tgj . The natural framework for discussing 
such corrections is chiral perturbation theory and we dis- 
cuss these chiral corrections next. 



II. CHIRAL PERTURBATION THEORY 
RESULTS FOR NUCLEON MAGNETIC 
MOMENTS AND MASSES 

In recent years there has been tremendous progress in 
the calculation of hadron properties using lattice QCD. 
Moore's Law, in combination with sophisticated algo- 
rithms, means that one can now make extremely accurate 
calculations for light quark masses (m^) larger than 50 
MeV. However, in order to compare with experimental 
data, it is still necessary to extrapolate quite a long way 
as a function of quark mass. This extrapolation is ren- 
dered non-trivial by the spontaneous breaking of chiral 
symmetry in QCD, which leads to Goldstone boson loops 
and, as a direct consequence, non-analytic behaviour as 
a function of quark mass [21I |2^ . Fortunately the most 
important nonanalytic contributions are model indepen- 
dent, providing a powerful constraint on the extrapola- 
tion procedure. 

In the past few years the behaviour of hadron proper- 
ties as a function of quark mass has been studied over 
a much wider range than one needs for the present pur- 
pose m m m in m Ei . one can therefore apply 
the successful extrapolation formulas developed in the 
context of lattice QCD with considerable confidence. 

The key qualitative feature learnt from the study of 
lattice data is that Goldstone boson loops are strongly 
suppressed once the Compton wavelength of the boson 
is smaller than the source. Inspection of lattice data for 
a range of observables, from masses to charge radii and 
magnetic moments, reveals that the relevant mass scale 
for this transition is rUq ~ 50 MeV - i.e., m^r ^ 400 — 500 
MeV [23. I29I I. The challenge of chiral extrapolation 
is therefore to incorporate the correct, model indepen- 
dent non-analytic behaviour dictated by chiral symme- 
try while ensuring excellent convergence properties of 
the chiral expansion in the large mass region, as well 
as maintaining the model independence of the results of 
the extrapolation. Considerable study of this problem 
has established that the use of a finite range regulator 
(FRR) fulfils all of these requirements js^, IS |33| ■ In- 
deed, in the case of the mass of the nucleon, it has been 
shown that the extrapolation from rri^ ^ 0.25 GeV2 to 
the physical pion mass - a change of by a factor of 
10 - can be carried out with a systematic error less than 
1% HJl- In the following we apply this same method to 
calculate the change in the nucleon mass, corresponding 
to quark mass changes at the level of 0.1% or less, as 
required in the present context. 



A. Variation of the nucleon mass with quark mass 

The expansion for the mass of the nucleon given in 
Ref. [3II33 is: 

Mm = ao-l-a2m^-f a4m^-f aem^+o-AT^+o-ATT+o-tad , (8) 
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where the chiral loops which given rise, respectively, to 
the leading and next-to-leading nonanalytic (LNA and 
NLNA) behaviour are: 
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32 7r/2 



g\ hi (tott , Aatat , A) (9) 



3 32 



Ctad 



32 7r/2 25 



gA^Af(m^, AjvA,A) (10) 



C2m^/T(m^, A) 



(11) 



and the relevant integrals are defined (in heavy baryon 
approximation) as: 



lM{mp,ABB',J^) = - / 
TT Jo 

/T(m7r, A) : 



dk^;^^^ (12) 



dk 

Xy^k'^ + m 



to , (13) 



with ujk — i/fc^ -|- mp and Abb' the relevant baryon 
mass difference (i.e., Mb' — A/s). We take the A~A^ 
mass splitting, A — Ma ~ M^f, to have its physical value 
(0.292 GeV), while qa = 1.26. The regulator function, 
u{k,A), is taken to be a dipole with mass A = 0.8 GeV. 
In Eq. H13|l to, defined such that It vanishes at rriTr = 0, 
is a local counter term introduced in FRR to ensure a 
linear relation for the renormalisation of C2. 

The model independence of the expansion given in 
Eq. ||HJ) is ensured by fitting the unknown coefficients 
to the physical nucleon mass and lattice data from the 
CP-PACS Collaboration jH, yielding: qq = 1.22, aa = 
1.76, a4 — — 0.829, ag — 0.260 (with all parameters ex- 
pressed in the appropriate powers of GeV). With these 
parameters fixed one can evaluate the rate of change of 
the mass of the nucleon with quark or pion mass at the 
physical pion mass: 



d 



drua 



Mn — m 



d 



■M 



N 



0.035 GeV, 



(14) 



a quantity commonly known as the pion-nucleon sigma 
commutator. Using Eq. (|14|l one finds the relationship 
(in terms of dimensionless quantities): 



6Mn m"^ OMn Sruq 
Mn dm^ ruq 



Mn 



0.037^ 



(15) 
(16) 



The extension of this procedure to the effect of a vari- 
ation in the strange quark mass is similar, but one must 
include the variation arising from 77- Nucleon loops, as 
well as Kaon loops with intermediate S or A baryons. 



'NT. 



'NN 



These contributions can be expressed as 
3 



p 

<^BB' 



327r/2 



Gbb' lM{'rnp,ABB',A) 



(17) 



(18) 



with G^g, the associated coupling squared. Once again 
we select the dipole regulator: 



u{k,A) = 



A2 



A2 + fc2 



For the relevant diagrams, N YiK , N 
N — > Nr], we have 



(19) 



KK and 



G 



G 



K 

NT, 

K 
NA 



'^NN 



{3F + D)^ 
{iF-Df 



(20) 



where we take F = 0.50 and D = 0.76. We use the Gell 
Mann-Oakes-Renner relation in the SU(2) chiral limit to 
relate the variation of the kaon mass in the chiral SU(2) 

limit, = 



= 0.484 GeV (with ^i^{K], the 



physical pionjkaon} mass) , to the variation of the strange 
quark mass {5in\/ih\ = Sms/iris). Hence the variation 
of the nucleon mass with strange quark mass is given by: 



N 



N 



d 



'NT 



^NA 



'NN/ 



(21) 

Using the dipole regulator mass, A = 0.8 GeV, Eq. H21() 
leads to the result: 



^ = 0.011^ 



(22) 



B. 



Variation of proton and neutron magnetic 
moments with quark mass 



The treatment of the mass dependence of the nu- 
cleon magnetic moments is very similar to that for the 
masses. Once again the loops which give rise to the 
LNA and NLNA behaviour are evaluated with a FRR, 
while the smooth, analytic variation with quark mass is 
parametrized by fitting relevant lattice data with a finite 
number of adjustable constants. 

For the lattice data we use CSSM Lattice Collabora- 
tion results 's^l of nucleon 3-point functions. Results are 
obtained using established techniques in the extraction 
of form factor data 35]. Similar calculations have also 
been recently reported by the QCDSF Collaboration [28^ . 
We use the two heaviest simulation results, ~ 0.6- 
0.7 GeV^ HI. These simulations were performed with 
the FLIC fermion action |3^ on a 20^ x 40 lattice at 
a = 0.128 fm. 

In the magnetic moment case the formulae are a little 
more complicated, so we leave the details for the Ap- 
pendix. Suffice it to say here that the relevant processes 
are shown in Fig. ^ Again we use a dipole form for the 
regulator with A = 0.8 GeV. 
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(a) 



(b) 



(d) 
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FIG. 1: Chiral corrections to the nucleon magnetic moments 
included in the present work. 



Having parametrized the neutron and proton magnetic 
moments as a function of m^, the fractional change ver- 
sus nriq or rUs is given by: 



5 PL ( d ^J■ \ SrUg 
fi \ dml } mq 



(23) 



A* 



(24) 



The numerical results may then be summarised as: 



^J'P 



-0.087^ 



Hp rus 



Pn 



Sm„ 
-0.118 2 



^ = 0.0013^ 



(Pp/fJ-n) rriq 



^(MpMO ^ _Q 015— 



{ppl fin/ 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



III. DEPENDENCE OF ATOMIC TRANSITION 
FREQUENCIES ON FUNDAMENTAL 
CONSTANTS 

Using the results of the previous section we can now 
use Eqs. (|6I7|) to study the variation of nuclear magnetic 
moments. For all even Z nuclei with valence neutron 



(i99jjg^i7iYb,iiiCd, etc.) we obtain ^ = ^ 
we have a valence proton with j—7/2, l^A and 

6m^ 



Sp 



0.110^ 



0.017- 



For Cs 



(31) 



For *^Rb we have valence proton with j=3/2, l—l and 



Sp 



-0.064^ 



0.010- 



(32) 



As an intermediate result it is convenient to present the 
dependence of the ratio of the hyperfinc constant, A, to 

4 

the atomic unit of energy E = (or the energy of the 
ls-2s transition in hydrogen, which is equal to 3/8 i?) on 
a variation of the fundamental constants. We introduce 
a parameter V defined by the relation 



SV_ ^ 5{A/E) 
'V ^ A/E 



(33) 



We start from the hyperfine structure of ^^^Cs which 
is used as a frequency standard. Using Eqs. H2I31|I we 
obtain 



,2.83 



A, 



QCD 



-jO.llO^ 



A 



QCD 



'■^ (34) 



The factor ^ will cancel out in the ratio of hyperfine 
transition frequencies. However, it will survive in com- 
parison between hyperfine and optical or molecular tran- 
sitions (see below). According to Eqs. (|16|) and l|22|) the 
relative variation of the electron to proton mass ratio can 
be described by the parameter 



X{me/mp) = ( 



A, 



QCD 



_-j-0.037j- 



A 



-) 



-0.011 



QCD 



A, 



QCD 



(35) 



which can be substituted into Eq. H34|l instead of nie/mp. 
This gives an expression which is convenient to use for 
comparison with optical and molecular vibrational or ro- 
tational transitions 



A 



QCD 



■|0.073j- 



\0.006 



A, 



QCD 



A, 



QCD 



(36) 



The dependence on the strange quark mass is relatively 
weak. Therefore, it may be convenient to assume that 
the relative variation of the strange quark mass is the 
same as the relative variation of the light quark masses 
(this assumption is motivated by the Higgs mechanism of 
mass generation) and to use an approximate expression 

F(133Cs) « a2-83(-^)0-13^. 

^ ' ^ Aqcd ' nip 

For hyperfine transition frequencies in other atoms we 
obtain 



2.34/ 



VC'Rb) = a 



m„ 



'A 



QCD 



_^-0.064|^ 



TO. 



A, 



-) 



-0.010 



m, 



QCD 



— (37) 



rrir, 



2( ™i? •j-0.087^_ 



A 



QCD 



■A. 



QCD 



-0.013 ™ 

m 



(38) 



p 



5 



VCH) = a" 



QCD 



_^-0.018|| 



A, 



QCD 



(39) following form 



^QCD ^QCD mp 



■^QCD ^QCD ' 



y("icd+) = a2-6(-:^)-°-"«(- 

(42) 

Note that the hyperfine frequencies of all even-Z atoms 
where the nuclear magnetic moment is determined by 
a valence neutron have the same dependence on quark 
masses. 



IV. LIMITS ON VARIATION OF 
FUNDAMENTAL CONSTANTS 

Now we can use these results to place limits on the pos- 
sible variation of the fundamental constants from partic- 
ular measurements. Let us start from the measurements 
of quasar absorption spectra. Comparison of the atomic 
H 21 cm (hyperfine) transition with molecular rotational 
transitions 9) gave limits for the variation of Yg = a^gp. 
In Refs. 5, 37] it was suggested that one might use these 
limits to estimate variation of mq/ Kqcd- According to 
Eqs. H25|) and (|26|) the relative variation of Yg can be 
replaced by the relative variation of Y {5Y/Y — 5Yg/Yg) 



Y 



"T-g ^ -0.087 1- ™s - 
^QCD ^QCD ' 



-0.013 



(43) 



Then the measurements in Ref. p lead to the following 
limits on the variation of Y: 

SY/Y = (-0.20 ± 0.44)10-5 for redshift z=0.2467 and 
SY/Y = (-0.16 ± 0.54)10^5 for z=0.6847. 
The second limit corresponds to roughly t=6 billion 
years ago. There is also a limit on variation of = 
a^gpiTie/mp obtained in Ref. |lOj |. This limit was inter- 
preted as a limit on variation of a or nif./mp. The relative 
variation of X„i can be replaced by the relative variation 
of 



i^QCD i^QCD A 



QCD 



(44) 



The dependence on quark masses appears from both the 
proton g-factor and the proton mass. The measurement 
in Ref. 10] leads to the following limit on variation of X: 
5X/X = (0.7 ± 1.1)10-5 for z=1.8. 

Now let us discuss the limits obtained from the labo- 
ratory measurements of the time dependence of hyper- 
fine structure intervals. The dependence of the ratio of 
frequencies A(-'^'^^Cs)/A(*^Rb) can be presented in the 



X{Cs/Rb) 



V{Cs) 



V{Rh) 

= a'-'^[mJKQCDr''[mslKQCDr^' (45) 

and the result of the measurement in Ref. may be 
presented as a limit on variation of the parameter X: 

1 dX{Cs/Rb) 
X{Cs/Rb) dt 



(0.2 ±7) X 10-'" /year. (46) 



Note that if the relation (QJ were correct, the varia- 
tion of X{Cs/Rb) would be dominated by variation of 
[niq/ Aqcd]- The relation (Q) would give X{Cs/Rb) oc 

For A(i33Cs)/A(H) we have 



X{Cs/H) 



V{Cs) 



V{H) 

= a^-^^[mJKQCDr''[ms/KQCDr'' (47) 

and the result of the measurements in Ref. may be 
presented as 

1 dX(Cs/H), 
^ XiCs/H) dt V 5.5xlO-^Vyea.. (48) 

For A(i99Hg)/yl(H) we have 

« a'-'[mg/AQCDr'-"''[ms/AQCDr'' ■ (49) 

The result of the measurement in Ref. [l5| may be pre- 
sented as 



1 



dX{Hg/H) 



'X{Hg/H) dt 



I < 8 X IQ-'^/year 



(50) 



Note that because the dependence on masses and strong 
interaction scale is very weak here, this experiment may 
be interpreted as a limit on the variation of a. 

In Ref. ^3 a limit was obtained on the variation of the 
ratio of hyperfine transition frequencies ^^^Yb"'"/^^'^Cs 
(this limit is based on the measurements of Ref. 38]). 
Using Eqs. (|34I41|I we can present the result as a limit 
on XiYb/Cs) = aO-7K/AQCD]-°-''*^K/AQCD]-°-°^': 



1 



dX{Yb/Cs) 



X{Yb/Cs) 



dt 



1(2) xlQ-'y year. (51) 



The optical clock transition energy E{Hg) (A=282 nm) 
in the Hg+ ion can be presented in the following form: 

4 

m e 

E{Hg) = const x [—^]FreiiZa) . (52) 

h 

Numerical calculation of the relative variation of E{Hg) 
has given 0: 



E{Hg) 



-3.2—. 

a 



(53) 
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This corresponds to V{HgOpt) = a~^-^ . Variation of 
the ratio of the Cs hyperfine spUtting A{Cs) to this 
optical transition energy is described by X{Opt) = 
V{Cs)/V{HgOpt): 



A 



QCD 



QCD 



'A 



QCD 



Here we used Eq. ^ for V{Cs). The work of Ref. |2C 
gives the hmit on variation of this parameter: 



dX{Opt) 



'X(Opt) dt 



< 7 X 10^^^ /year . 



(55) 



Molecular vibrational transitions frequencies are pro- 
portional to {me/rripY/'^ . Based on Eq. (|H5|I we may 
describe the relative variation of vibrational frequencies 
by the parameter 

V{vib) = (^!^)-0.018(^^)-0.005(^!^)0.5 (56) 

A.QCD Aocn At 



^QCD 



^QCD 



Comparison of the Cs hyperfine standard with 
SFq molecular vibration frequencies was discussed 
in Ref. 39]. In this case X{Cs /Vibrations) = 

0.091/ 



a'^-^[me/AQCD]°-^[mq/AQCD 



0.011 



The measurements of hyperfine constant ratios in dif- 
ferent isotopes of the same atom depends on the ra- 
tio of magnetic moments and is therefore sensitive to 
rriq/ Aqcd- For example, it would be interesting to mea- 
sure the rate of change for hydrogen/deuterium ratio 
where X{H/D) = K/AQCD]-*'-°^«[m,/AQCD]°-°3'. 

Walsworth has suggested that one might measure the 
ratio of the Zeeman transition frequencies in noble gases 
in order to explore the time dependence of the ra- 
tio of nuclear magnetic moments. Consider, for exam- 
ple ^^^Xe/'^He. For '^He the magnetic moment is very 
close to that of neutron. For other noble gases the nu- 
clear magnetic moment is also given by the valence neu- 
tron, however, there are significant many-body correc- 
tions. For ^^^Xc the valence neutron is in an si/2 state, 
which corresponds to the single-particle value of the nu- 
clear magnetic moment, [i = [in = —1.913. The mea- 
sured value is /i = —0.778. The magnetic moment of 
the nucleus changes most efficiently through the spin- 
spin interaction, because the valence neutron transfers 
a part of its spin, < Sz >, to the core protons and 
the proton magnetic moment is large and has the op- 
posite sign. In this approximation ^ — [1 — b)fj,n + 
bnp. This gives b=0.24 and the ratio of magnetic mo- 
ments y = Ai£f^Xe)/Ai(3He)« 0.76 -I- 0.24gp/g„. Using 
Eqs. (|25I26I27I28|I we obtain an estimate for the relative 
variation of /i(^^^Xe)//i(^He), which can be presented as 
variation of X = K/AQci3]~° *"^[m,/AQCD]°-"''. Here 
again 6Y/Y = SX/X. 

Note that the accuracy of the results presented in this 
paper depends strongly on the fundamental constant un- 
der study. The accuracy for the dependence on a is a few 
percent. The accuracy for ruq/AqcD is about 30% - be- 
ing limited mainly by the accuracy of the single-particle 



TABLE I: Chiral coefficients for various diagrams contribut- 
ing to proton and neutron magnetic moments. We use SU{Q) 
symmetry to relate the meson couplings to the ttA'^A vertex, 
C = -2D. 



a 




B" 


(a) 


-{F + Df 


{F + Df 


(b) 




9*- 


(c) 







(d) 


-\{D-Ff 


-{D-Ff 



approximation for nuclear magnetic moments. (For com- 
parison, the estimated systematic error associated with 
the calculation of the effect of the quark mass variation 
is less than 10%.) Finally, we stress that the relation ^ 
between the variation of a and m/ Aqcd has been used 
solely for purposes of illustration. 
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APPENDIX A: MAGNETIC MOMENTS 

As explained in the text, we explicitly include the pro- 
cesses shown in Fig. ^ which give rise to the leading and 
next-to-leading nonanalytic behaviour as a function of 
quark mass. 

We describe the quark mass dependence of the mag- 
netic moments as: 



a2ml 



(Al) 



where denotes the chiral loop corrections given by 



+XM(b) ^mC^tt, Aata, A) 

+Xm(c) Ajva, A) 

+X^,{d)If^{mK,^N^,A)- 



(A2) 



The chiral coefficients of the loop integrals, X/^a, are given 

by 



M 



N 



(A3) 



and are summarised in TableHlliiliJ,!!!. Note that the 
required analytic terms in the chiral expansion to this 
order have been placed in a Fade approximant designed 
to reproduce the Dirac moment behaviour of the nucleoli 
at moderate quark mass. 
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The corresponding loop integral is given by 



4 /■= 

/,KA,A) = --y^ 



dk 



(A + 2cjfe)fc4u2(fc,A) 



(A4) 



where the various terms have been defined in Sect. II. We 
note that in the limit where the mass-splitting vanishes 
this integral is normalised such that the leading nonana- 
lytic contribution is m. 

With the coeflicients of the loop integrals defined, we 
only require determination of the parameters ao and a2 
in Eq. ljAl|l to constrain the variation with quark mass. 
We note also that this form assumes no analytic depen- 
dence on the strange quark mass, beyond what is implic- 
itly included in the loop diagrams (c, d). We determine 
cto,2 for both the proton and neutron by fitting the phys- 
ical magnetic moment as well as the lattice QCD data. 
We fit only to the two heaviest simulation results of the 
CSSM Lattice Collaboration [H, ml ~ 0.6-0.7 GeV^ 
These simulations were performed with the FLIC fermion 
action (36j| on a 20^ x 40 lattice at a = 0.128 fm. We se- 
lect the heaviest two data points, where the effects of 
quenching are anticipated to be small |43l li^ . 



The best fits to the physical values and the lattice data 
give 

= 2.17 fj.N al = 0.817 GeV"^ (A5) 
= -1.33 "2 = 0.758 GeV"^ . (A6) 

Upon renormalisation of the loop diagrams, the resultant 
magnetic moments in the SU(2) chiral limit are given by 



/iQ — 3.48 , and /Iq — —2.58 /ijv ■ 



(A7) 



We now take derivatives of Eq. (|A1|) at the physical 
pion mass to determine the variation with quark mass. 
In particular, we have 



r ml d/i 1 5mq 



/i dm^ J TO, 



d II 1 Sm 



/i drh^ 



K 



(A8) 



(A9) 



This yields the results shown in the text. 
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